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Abstract. Analytical and numerical vortex solutions for the extended Skyrme-Faddeev model 
| in a (3 + 1) dimensional Minkowski space-time are investigated. The extension is obtained by 

adding to the Lagrangian a quartic term, which is the square of the kinetic term, and a potential 
which breaks the SO(3) symmetry down to SO{2). The construction of the solutions has been 
done in twofold: one makes use of an axially symmetric ansatz and solves the resulting ODE 
by an analytical and a numerical way. The analytical vortices are obtained for special form 
^ ' of the potentials, and the numerical ones are computed using the successive over relaxation 

£Nj , method for wider choice of the potentials. Another is based on a simulational technique named 

■ the simulated annealing method which is available to treat the non-axisymmetric shape of 

solutions. The crucial thing for determining the structure of vortices is the type of the potential. 

OV 
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1. Introduction 

The so-called Skyrme-Faddeev model was introduced in the seventies pQ as a generalization to 
(3 + 1) dimensions of the 0(3) non-linear sigma model in (2 + 1) dimensions. The Skyrme term, 
quartic in derivatives of the field, balances the quadratic kinetic term and according to Derrick's 
theorem, allows the existence of stable solutions with non-trivial Hopf topological charges. Due 
to the highly non- linear character of the model and the lack of symmetries, the first soliton 
solutions were only constructed in the late nineties using numerical methods [21 [31 SI E] . Since 
then the interest in the model has increased considerably and it has found applications in many 
areas of physics due mainly to the knotted character of the solutions [6j. One of the aspects of 
the model that has attracted considerable attention has been its connection with gauge theories. 
Faddeev and Niemi have conjectured that it might describe the low energy limit of the pure 
SU{2) Yang-Mills theory [7]. They based their argument on a decomposition of the physical 
degrees of freedom of the SU{2) connection, proposed in the eighties by Cho [8], and involving 
a triplet of scalar fields n taking values on the sphere S 2 (n 2 = 1). Gies [9] has calculated 
the Wilsonian one loop effective action for the pure SU{2) Yang-Mills theory assuming Cho's 
decomposition, and found that the Skyrme-Faddeev action is indeed part of it, but additional 
quartic terms in the derivatives of the triplet n are unavoidable. In fact, the first numerical 
Hopf solitons were first constructed for the Skyrme-Faddeev model modified by a quartic term 



[2] which is the square of the kinetic term. However, the soliton solutions in [2] were constructed 
for a sector of the theory where the signs of the coupling constants disagree with those indicated 
by Gies' calculations. Therefore, it is worth investigating the model with correct sign of the 
coupling constants. 

In this paper we consider an extended Skyrme-Faddeev model (ESF) defined by the 
Lagrangian 

C = M 2 duH ■ d^n - \(dnfi A d v nf + - (dji ■ d*n) 2 - V(n 3 ) (1) 
e / 2 

where n is a triplet of real scalar fields taking values on the sphere S 2 , n 3 its third component, 
M is a coupling constant with dimension of (length) -1 , e 2 and /3 are dimensionless coupling 
constants, and the potential is a functional of the third component n 3 of the triplet ft. Note 
that the potential breaks the 0(3) symmetry of the original Skyrme-Faddeev down to 0(2), 
the group of rotations on the plane n\ ni, and so eliminating two of the three Goldstone boson 
degrees of freedom. In this paper the main role of potential is to stabilize the vortex solutions. 

The static energy density (^static = — £) associated to (HJ is positive definite if V > 0, 
M 2 > 0, e 2 > and /3 < 0. That is the sector explored in [2] and where Hopf soliton solutions 
were first constructed (for V = 0). In addition, that is also the sector explored in [TU] but with 
additional terms involving second derivatives of the n field, and where Hopf soliton were also 
constructed. The static energy density of ([T]) is also positive definite for V > if 

M 2 > ; e 2 < ; (3 < ; /3 e 2 > 1 (2) 

That is the sector that agrees with the signature of the terms in the one loop effective action 
calculated in [Qj and it is the sector that we will consider in this paper. Static Hopf solitons 
were constructed in [11] for the sector ([2]) (with V = 0) and their quantum excitations, including 
comparison with glueball spectrum, were considered in |12j . An interesting feature of the Hopf 
solitons constructed in [11] is that they shrink in size and then disappear as (3 e 2 — > 1, which is 
exactly the point where the exact vortex solution exist [13] . 

The aim of the present paper is to investigate if vortex solutions for the model ([1]) continue 
to exist when the condition j3 e 2 = 1 is relaxed, and so if they co-exist with the Hopf solitons of 
|llj . In order to stabilize the solution, we shall introduce the types of potential 

V(n 3 ) = ^vl v ^(l + n 3 y(l-n 3 ) b (3) 

where a + b = non-zero integer, and \x is a real coupling constant. A special choice of the 
parameters a, b we have holomorphic solutions of the model while for the other case we still have 
numerical solutions. 

In this paper, first we discuss the integrable holomorphic solutions of the model and next we 
shall perform the numerical stuff. The numerical simulations are done by twofold: one is by 
solving a differential equation which is accomplished in terms of the standard successive over 
relaxation. Another is based on energy minimization scheme called the simulated annealing. 
Especially the latter analysis demonstrates the detailed behavior of the symmetry breaking of 
the solutions by the change of the structure of the potential. 

2. The integrable sector of the model 

The first exact vortex solutions for the theory ([1]) were constructed in [13] for the case where the 
potential vanishes, and by exploring the integrability properties of a submodel of ([T]). In order 
to describe those exact vortex solutions it is better to perform the stereographic projection of 



the target space S 2 onto the plane parameterized by the complex scalar field u and related to n 
by 



n = (u + u*, -i(u - u*),u 2 - 1)/(1 + u 2 ) (4) 
It was shown in [13] that the field configurations of the form 

u = u(z,y), u* = u*(z*,y), for f3e 2 = 1, V = (5) 

are exact solutions of ([T]), where z = x 1 + ieix 2 and y = x ?J — E2X , with e a = ±1, a = 1,2, 
and 2;^, /i = 0, 1,2,3, are the Cartesian coordinates of the Minkowski space-time. The simplest 
solution is of the form u = z n e lky , with n integer, and it corresponds to a vortex parallel to the 
x 3 -axis and with waves traveling along it with the speed of light. 

In terms of the complex scalar field u introduced in H] the Lagrangian Q] becomes 



, „ r2 d u ud^u* 8 

4M — - 1 

(1 + u 2 ) 2 ^ e 2 



(c» 2 (cU*) 2 2 Ad.ud^u*) 2 
(l + « 2 ) 4 {l ' {l + u 2 Y 



V(\ u | 2 ) (6) 



where we have used the fact that 713 is a functional of | u \ 2 only, and so is the potential. The 
Euler-Lagrange equations following from ([6]), or dU), reads 



(1 + M 2 )d M /C M - 2u*K^u = + |-u| 2 ) 3 V (7) 



where V = and 



4 \{d v ud v u)dau* + {l3e 2 -l)(d v ud u u*)dau] 

K„ = M 2 d^u + - ^ Lt W ———— ULl 8 

e z (1 + u z y 



We point out that the theory © possesses an integrable sector defined by the condition 

(d^u) 2 = (9) 

Such condition was first discovered in the context of the CP 1 model using the generalized zero 
curvature condition for integrable theories in any dimension [14] , and then applied to many 
models with target space being the sphere S 2 , or CP 1 (see |15| for a review). It leads to an 
infinite number of local conserved currents. Indeed, ([9]) together with the equations of motion 
([7]) imply the conservation of the infinity of currents given by 

where G is any functional of | u | 2 only. For the case where the potential vanishes, the set of 
conserved currents is considerably enlarged since G can be an arbitrary functional of u and u*, 
but not of their derivatives. If in addition to the condition Q one takes V = and f3 e 2 = 1, 
then the equations of motion reduce to d 2 u = 0. It is in that integrable sector that the solutions 
© lie, and were studied in [13]. 

It is interesting to note that ([7]) with a special choice of the potential 

y(n 3 ) = ^(l + n3) 2 - 2 / n (l-n 3 ) 2 + 2 /" (11) 
have an analytical, holomorphic solution for each topological charge as 

U(p, if, Z, T) = (£) n e il™<P+k(z+r)] ( 12 ) 



Figure 1. The n = 2 profile g(y) and the corresponding Hamiltonian density of the real space 
H(p) of k 2 = 0.0 for the constant r 2 p 2 /M 2 = 1.0. 



where e = ±1 and a describes a scale of the solution. Here we used dimensionless polar 
coordinates (p, (p, z, r) defined by 



x° = ct = tqt, x 1 = tqp cos(/9, x 2 = r$p siny?, x 3 

and where we have introduced a length scale r$ given by 

4 
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Substituting (fl2j) into (|7|) the a can be determined such as 



a = n 



M 2 (f3e 2 



rip 2 



1/4 



n 



(-e 2 ) (/3e 2 - 1)M 4 
4p 2 



1/4 



(15) 



Clearly, the special solution at fie 2 = 1 is obtained if we take a proper limit of the vanishing 
potential, i.e./3e 2 — > 1 and p 2 — >■ with (f3e 2 — l)/p 2 = constant. 

In the case with the potential, the current (|10p is still conserved because 



(16) 
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where we have used the reduced integrable equation 

• 2 ,dv{\u\ 2 ^ 



d^u*Kl = 



+ ^-(1 + \u 



|2\2 ' 



9u* 



0. 



(17) 



The Hamiltonian density associated to © is not positive definite due to the quartic terms in 
time derivatives. We shall arrange the Legendre transform of each term in §6§ to make explicit 





Figure 2. The static energy and its components corresponding to the solutions of Fig{TJ 



such non positive contributions, and write the Hamiltonian density as (see [16] for details) 



U 



\M 



u I 2 +Vu ■ Vu* 
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24 (Vu) 2 (Vu*) 



- 24 
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+ V(\u\ 2 ) 



(18) 



:i + n 2 ) 4 

where u denotes the x°-derivative of u, and Vu its spatial gradient, and where we have denoted 

1 



u 2 



(Vu)< 



+ Fe l 



(19) 



with F > and < < 2tt, being functions of the space-time coordinates. The most of terms 
in (|18p make positive contribution while the second term has some possibility to be negative. 
Note also, for static configurations apparently it is positive definite for the range of parameters 
given in ([2]). 

3. The integrable and the non-integable sectors: numerical analysis by the SOR 

Although in the previous section we used the polar coordinates, for the numerical study it is 

more convenient to use a new radial coordinate y (0 < y < 1), defined by p = \J^~^~- We 
introduce the solution ansatz of the form 



(20) 



where the profile function g{y) is defined at the period < g < 1. The equation can be written 
as 



d_ 

dy 
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where the primes at this time indicate derivatives w.r.t.y and where 

n = n 2 



(21) 



S = l + Pe 2 g{l- g )^ + (u^ g i 
R = 1 + g{l - g)^L | (/3e 2 - 2)0 + /3e 2 (j^g* 



(22) 



The energy in the unit of 4M 2 per unit length for the time-dependent vortex can be estimated 
in terms of following four parts of integrals of the dimensionless Hamiltonian H := 71 /AM 2 



roc 2 2 

E = 2tt pdpH(p) =E 2 + Efp + ((3e 2 - l)i?f + 7 j^E [ 



(23) 




Figure 3. The n = 2 static energy 
k 2 = and r^ 2 /M 2 = 1.0. 
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and the component E± for several type of potentials 



v ■ 



in which the components are defined as 
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(26) 
(27) 



It is easy to see that the k 2 term in (|26p has positive contributions to the energy while in (|25p 
the sign of the fe 2 term depends on the spatial structure of the solution g. For the holomorphic 
solution (fl2j) . the energy ([25]) becomes zero and then the energy of the integrable sector keeps 
positive definite for all values of k 2 . 

For n = 2, the explicit form the potential is 



V, 



n=2 



;i + n 3 )(l-n 3 ) 3 . 



(28) 



The potential has zero at both the origin and the infinity thus it is so called a new-US (Baby- 
Skyrmion) type. The equation ([2~T]) . (f22|) are solved in terms of the standard successive over 
relaxation scheme. Fig(TJ is the profile function and the Hamiltonian density for n = 2. Figj2] 
is the energy per unit length and its components for several values of fie 2 and fixed /j, 2 . We 
confirmed that the value of the component e[^ is regarded as zero within the numerical 
uncertainty. This clearly indicates that the solution satisfies the condition ([9]). 

Although we have obtained the analytical solutions for a special form of the potential ([3]), 
we have many options for choice of the potential. We can obtain many numerical solutions for 



Figure 4. The energy density per unit length of the oW-BS potential v® (left) and the new-BS 
potential v\ (right). 




the several types of the potentials. We show the result of n = 2 for the potentials w>, v®, v g ^ 3 ; 
of course these are not of the form of the analytical solution. Note that for such non-integrable 
solutions, the second term in (|18p survives and has a possibility to be negative. Here we only 
consider the case for k 2 = 0. FigJ3] presents the energies and the component E± for these 
potentials. For n = 2, the old-US potentials give higher total energy than the new-BS. This 
indicates that the same class of potentials gives the similar energy and then, for n = 2 the 

4/3 

energy of the new type potential t> g ^ 3 is closest to the integrable sector, which is also plotted in 
Fig C2 for reference. 



4. Broken axisymmetric solutions: analysis by the SA 

In the previous sections, we have assumed that the solution is invariant under the 0(2) internal 
symmetry u — > e^u as well as the transformations of the Poincare group given by rotation on 
the plane x l x 2 and translations in the directions x° and x 3 . However, there is a possibility of 
relaxing some of these symmetries. In [T7], the authors generalized the old-BS type potential as 

2 

V = ^-(1 — nsY (0 < s < 4) and saw how the solution deforms from the rotational symmetry 
on the x 1 x 2 plane by the change of fi or s. 

For the deformed solution, the straightforward generalization of the ansatz ()20p is 



u = f(p,ip)e ieM e ikiz+r) 



(29) 



Figure 6. The energy density per unit length of V e = ^-(1 — w 3 )(l + era 3 ) for e = 0.25 (left) 
and e = 0.5 (right). The results for e = 0, 0, 1.0 are already presented in FigHl 

or equivalently, 

n = (sinF(/), p) cos[G(/9, ip) + k(z + r)], sinF(p, if) sin[©(/3, p) + A;(z + r)], cos -F(/9, <^)) (30) 

where the n field in terms of u is given in Q. Here ©(p, </>) := rap + ©o(p, V 3 ) an d ©0 is 
homotopic to the constant map. The method is a kind of the Monte-Carlo method in which 
one generates random numbers and properly change the value of the fields n by the numbers so 
as to drop the energy. However, a more sophisticated method may be applied to the problem. 
The simulated annealing method [18J is the application of the Metropolis algorithm which can 
successfully avoids the unwanted saddle points. The analysis has been done by minimizing the 
energy per unit length which is obtained by integrating the hamiltonian corresponding to (UJ 
with the ansatz (|30p into the (p, p) plane. We have confirmed that for the holomorphic solution 
(|12p the hamiltonian is positive definite. For the general solution, however, the energy positivity 
is supported only for a certain range of k. Here we present the results for k = 0. Also, we 
examine the case of n = 3 because the solution exhibits obvious symmetry breaking such as 
from axial to Z2-symmetry. 

We study the following several cases: 

(i) We start with two standard cases: the old-BS potential v J and the new^BS potential v\. In 
FigJH we present the energy density per unit length. For the oW-BS potential, the solution 
strongly deforms from the axial symmetry while for the new-BS the solution keeps the 
symmetry. 

(ii) Next we shall see how the solution behaves for the change of the model parameters. In 
Fig03 we plot the energy density per unit length of the old-BS potential for several values of 
the model parameter j3e 2 . As is easily observed that for larger value of (3e 2 , the deformation 
is enhanced. 

(iii) From the results of (i),(ii), we get a new insight for the symmetry breaking. That is, most 
crucial thing for the deformation is that the potential is finite or not at the point antipodal 
to the vacuum. In order to investigate this criterion further, we introduce following two 
types of one-parameter family of potentials 

u 2 

K(n 3 ) := y(l-n 3 )(l + en 3 ), 0<e<l (31) 
Vs(n 3 ) := ^(l-n 3 )(l+n 3 ) 5 , 0<S<1 (32) 




Figure 7. The energy density per unit length of Vg = ^-(1 — ns)(l + n^Y for 5 = 0.25 (left) 
and <5 = 0.5 (right). The results for 5 = 0.0, 1.0 are already presented in FigHl 

Note that for e,8 = the potentials become the oW-BS and e,5 = 1 they are the netoBS. 
The results are stimulating. In the case of V e , the solutions always break the axial symmetry 
except only for e = 1.0 (FiglH]). On the other hand, in the V$, the solutions always keep 
the axial symmetry except only for 5 = 0.0 (FigJTJ). Thus we confirmed the criterion for 
the mechanism of the symmetry breaking: if a potential has the vacuum value at the point 
antipodal to the true vacuum, the solution always exhibits the axial symmetry, and if there 
is no another zero except for the vacuum, the solution deforms. 

We summarize our results: for n = 1, the potential for the integrable sector is the old-US type 
and the solution naturally has the axial symmetry. For n > 2 the potential for the integrable 
sector is always the new-US type, and in terms of the above simulational study, the solutions 
always should be axially symmetric. As a result, our prescription of the ansatz (|20p is valid for 
all topological charges. 
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